Skyrmions in the presence of isospin chemical potential 
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Abstract 

We analyze the existence of localized finite energy topological excitations on top of the pertur- 
bative pion vacuum within the Skyrme model at finite isospin chemical potential and finite pion 
mass. We show that there is a critical isospin chemical potential fi'f above which such solutions 
cease to exist. We find that fij is closely related to the value of the pion mass. In particular for 
vanishing pion mass we obtain fij = in contradiction with some results recently reported in the 
literature. We also find that below fj,j the skyrmion mass and baryon radius show, at least for the 
case of the hedgehog ansatz, only a mild dependence on the isospin chemical potential. 

PACS numbers: 12.39.Dc, 25.75.Nq 
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Hadronic systems with vanishing baryon chemical potential /x b and finite isospin chemical 
potential /i/ are unstable with respect to weak decays. However, if we are interested in the 
dynamics of the strong interaction alone, we can disregard the relative slow electroweak 
effects and consider them as stable. Moreover, although there are not yet precise lattice 
QCD calculations at finite baryon density due to the Fermion sign problem, it is in principle 
possible to perform lattice simulations at finite isospin density[l|. These remarks have led 
several groups to study the behavior of strongly interacting matter at /j,b — and finite /i/. 
Effective lagrangian analysis showed that there is a phase transition from normal phase to 
pion superfluidity at a critical isospin chemical potential which turns out to be equal to the 
pion mass in the vacuuml2|. This has been confirmed by lattice QCD calculations [3|, random 
matrix method analysis jj], etc. Studies at finite temperature have been also performed |5|. 
Given these results it is of considerable interest to investigate on the behavior of baryon 
properties at finite isospin chemical potential. One model which is well suited to perform 
these studies is the Skyrme model. In the Skyrme model 6( and its generalizations, baryons 
arise as topological excitations of a non-linear chiral Lagrangian written in terms of meson 
fields. These type of models have been quite successful in describing the properties of octet 
and decouplet baryons (see e.g., Refs.pJ). In a series of recent articles^ the skyrmion 
properties in the presence of the isospin chemical potential have been analyzed. It has 
been found that, in the case of vanishing pion mass, there is a critical chemical potential 
lij m 223 MeV above which stable soliton solutions cease to exist. Moreover, according to 
Refs. [8 1 the skyrmion mass vanishes at fi = /jj. In this letter we re-examine the issue of the 
skyrmion stability for finite isospin chemical potential considering the possibility of having a 
finite pion mass. We show that although there is indeed a critical isospin chemical potential 
/if, its value is closely related to the value of the pion mass, and for vanishing pion mass 
one has fj,j = 0. Moreover we find that, at least within the spherically symmetric hedgehog 
ansatz, the skyrmion mass and baryon radius remain rather stable up to the critical isospin 
chemical potential. 

We start by considering the lagrangian of the 577(2) Skyrme model with quartic term 
stabilization and finite pion mass. It is given by 

C = -^Tr {L a L a } + — Tr { [L a , L,f) + -^Tr {U + C/+ - 2} . (1) 

Here, f n is the pion decay constant whose empirical value is f% mp = 93 MeV, e is the 
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so-called Skyrme parameter and m n is the pion mass which we will take at its empirical 
value m^ mp = 139 MeV. In our numerical calculations below we will use the standard set 
of values f n = 54 MeV and e = 4.84 which leads to the empirical value of the nucleon and 
A masses within the rigid rotor approximation [9]. It should be stressed, however, that our 
main conclusions are expected to be independent of this particular choice of parameters. In 
Eq.(0Q), as usual, U represents the 577(2) chiral field and the Maurier-Cartan operator L a is 
defined by L a = Wd a U. 

The isospin chemical potential \x\ is introduced by performing the replacement 

d a U -^d a U-i^ [r 3 , U] g a0 , (2) 

where g a p is the metric tensor in Minkowski space and r 3 is the third Pauli matrix. This 
leads to a modified lagrangian which reads 

CM =C + ^Tr {u 2 } - ^Tr { [u, L a } 2 } + ^Tr {uL } - ^Tr {uL a [L , L a }} (3) 

where u = UW 3 U — r 3 . 

In what follows we will be interested in static soliton configurations. The corresponding 
soliton mass reads 

c2 1 ^2 £2 
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It is not hard to see that the terms proportional to fij are not invariant under isospin 
rotations. Namely, the isospin chemical potential introduces a preferred direction in isospin 
space which is expected to lead to an axially deformed soliton configuration. For the time 
being we will assume that such deformations are small for the range of values of fii considered 
here. Therefore we introduce the usual spherically symmetric hedgehog ansatz for the baryon 
number B = 1 configuration 

Uh = exp [if ■ f F(r)] . (5) 

In this case we obtain 

f 2 r r s 2 ( f' 2 \ 1 s 4 

MhM = fj^ \F' 2 + 2 72 [l + -TJl) + ^ +2ml(l-c) 



where s = sinF and c = cosF. The minimization of Mh(/zj) leads to the following Euler- 
Lagrange equation for the soliton profile F 
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As usual this differential equation is supplemented by the boundary conditions corresponding 
to B = 1 topological excitations on top of the perturbative pion vacuum, F(0) = tc and 
-F(oo) = 0. The associated baryon radius is given by 

/ O roo \ 1/2 

r B M=(— drr 2 s 2 F') . (8) 
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Before presenting the numerical soliton solutions we will analyze the behavior of the 
profile F(r) for large distances. Given the boundary conditions we can linearize Eq.(j7]) in 
that limit. We obtain 

r+?£-(»i-lri + l)F = 0. (9) 
This equation implies that localized finite energy topological excitations on top of the per- 



turbative pion vacuum exist only for /ij < where /if = a/3/2 m n . For fij > jjj the 
solutions have an oscillatory behavior. This is clearly displayed in Fig.l where we plot the 
numerical solutions of Eq.fjTJ) for some chosen values of In fact, the situation is quite 
similar to the one found long time ago in the study of the stability of the skyrmion under 
spin-isospin rotations [lol]. Note that for /// = the solution displays a 1/r 2 large distance 
behavior which is typical of the localized pion massless case. 

In the region where localized solutions exist we can study the behavior of the soliton mass 
and baryon radius as a function of the isospin chemical potential. This is shown in Fig.2. We 
observe that as /// increases, the soliton mass decreases while the radius increases. However, 
these effects are not too large (less than 5 % at \ij = fij), and Mff(/zj) never vanishes in 
such region. 

As mentioned at the beginning of this letter, general arguments indicate that in the me- 
son sector there is a phase transition from the normal phase (perturbative pion vacuum) to 
a pion condensed phase at jjq = m n . Comparing with the result obtained above we note that 
in the soliton sector it appears an extra factor a/3/2 in the corresponding critical value. A 



similar factor has been found in the study of the stability of the hedgehog skyrmion under 



spin-isospin rotation 



11] . As in that case, it reflects the fact that in the hedgehog approx- 



imation possible pion excitations are assumed to be spherical while, as already mentioned, 
the presence of the isospin chemical potential is expected to induce axially symmetric de- 
formations. In order to account for this fact we introduce a general axial ansatz which, in 
cylindrical polar coordinates (p, </>,z), is given byjlll. Il2|. 
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ip 3 + iT 3 i/j 2 + iipi (tl cos (f) + t 2 sin < 



(10) 



Here, ip a are the components of a unit vector ip(p,z) that is independent of the angular 
variable 0. The boundary conditions for finite energy solutions are that ip — > (0, 0, 1) as 
p 2 +z 2 — * oo and that on the symmetry axis p = the equations ipx = and d p ip2 = d p ip3 = 
must be satisfied. Using this ansatz, the mass for the static soliton configuration reads 
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where i = p, z. Minimizing this expression we obtain a set of coupled equations for the 
two independent functions that we take to be tpi(p,z) and ip2{p,z). Unfortunately the 
resolution of these equations implies a rather time-consuming numerical task. From the 



results obtained in the case of spin-isospin rot at ions [ll|, soliton properties are not expected 
to be too much affected by the axial deformations provided pi < pj. Thus we postpone this 
numerical analysis for a future work. It is important, however, to consider the linearized 
form of the equations for ipi and ip2, which are valid at large distances. They read 
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Thus, once axially symmetric configurations are considered, we have indeed p\ = m v . 
Namely, for pi < p c j the behavior of the deformed skyrmion mass and radius as a function of 
the isospin chemical potential is expected to be very similar to that shown in Fig.2, the only 
important difference being that the curve will end at pi/m^ = 1 instead of pi/m^ = a/3/2. 

In conclusion, we have re-examined the behavior of the skyrmion properties as a function 
of the isospin chemical potential pi. We have found that there is, indeed, a critical value 



of hi above which localized finite energy topological excitations on top of the perturbative 
pion vacuum cease to exist. Such critical value is closely related to the value of the pion 
mass. For the spherically symmetric hedgehog ansatz we find fxj = a/3/2 m^, while when 
axially symmetric deformations are allowed we obtain /if = as expected from general 
arguments in the meson sector. These results disagree with previous works[8] which report 
fi c j « 223MeV in the chiral limit (m n = 0). Moreover, for fij < fi c j we find only a mild 
dependence of the soliton mass and radius on fij. In particular, in contrast to the results 
in Ref-j^l our values of the soliton mass do not vanish for any value of /// < [l c j. It should 
be noted that these analyses of the soliton mass and radius dependence on /// have been 
performed using the spherically symmetric hedgehog ansatz. Although for that range of 
values of /i/ the axially symmetric deformations are expected to have only a minor effect 
on this dependence, the numerical resolution of the corresponding equations is required in 
order to confirm this expectation. This will also allow us to explore the region fij > We 
hope to report on these issues in forthcoming publications. 
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FIG. 1: (Color online) Hedgehog profiles F(r) for various values of the chemical potential hi. Note 
that localized solutions only exist for ni/nj < 1. 
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FIG. 2: (Color online) Soliton mass (full line) and baryon radius (dashed line) as a function 
of the isospin chemical potential [ii for the spherically symmetric hedgehog ansatz and taken 
with respect to their corresponding values at fij = 0. For our parameter choice such values are 
M H (0) = 864 MeV and r B (0) = 0.68 /m[9|. 
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